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SUMMARY  AND  CONCLUSIONS 

This  report  analyzes  the  effects  of  antenna  elevation  and  inclina- 
tion on  the  structure  of  air-to-air  VLF/LF  communication  signals  for 
frequencies  between  15  and  40  kHz.  Considered  are  1)  propagation 
broadside  to  the  horizontal  component  of  trailing-wire  transmitting 

antennas,  and  2)  ambient  and  moderately  disturbed  daytime  ionospheric 

-3  -5 

conditions  and  ground  conductivities  of  10  and  10  mhos/m.  The 
calculations  use  full-wave  methods,  and  fully  account  for  the  vertical 
inhomogeneity  of  the  ionosphere  and  the  curvature  of  the  earth. 
Numerical  results  are  given  for  the  excitation  factors,  height-gain 
factors,  and  attenuation  rates  of  the  first  few  TM  and  TE  modes.  Mode 
summations  are  carried  out,  and  detailed  plots  are  given  for  each  field 
component  versus  distance  for  antenna  altitudes  between  0 and  40  kft. 

For  ground-based  antennas,  the  analysis  shows  that  the  vertical 
electric-field  component  associated  with  TM  modes  is  dominant,  even  for 
electric-dipole  antennas  inclined  as  little  as  10°  with  respect  to  the 
horizontal.  For  a moderate  spread-debris  nuclear  environment,  the 
vertically  polarized  component  dominates  the  horizontally  polarized 
component  for  all  reasonable  antenna  inclinations,  antenna  elevations, 
and  ranges  beyond  about  1500  km.  This  dominance  occurs  because  of  the 
high  attenuation  suffered  by  TE  modes  when  ionospheric  reflection- 
heights  are  severely  depressed. 

For  ambient  conditions — and  presumably  for  weakly  disturbed 
conditions — the  horizontally  polarized  components  associated  with  TE- 
mode  propagation  are  dominant  for  antenna  elevations  of  25  kft  or  more, 
and  antenna  inclinations  smaller  than  about  20°,  measured  from  the 
horizontal. 

1 


A 


I 


CONTENTS 


SUMMARY  AND  CONCLUSIONS  1 

PREFACE  2 

SECTION  I.  INTRODUCTION  6 

II.  ELECTRIC  DIPOLE  FIELDS  7 

Excitation  Factors  8 

Height-Gain  Factors  12 

Attenuation  Rates  18 

Field  Strength  versus  Distance  21 

III.  DISCUSSION  35 

APPENDIX:  MATHEMATICAL  SUMMARY  36 

REFERENCES  Aft 


f 


FIGURE  1.  Excitation  factor  vs.  frequency  for  lowest  three 


TM  modes  (a-1,2,3)  and  lowest  two  TE  inodes 

(8-1,2):  ambient  day,  - 10“  3 mhos/m  10 
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I.  INTRODUCTION 


Propagation  of  very-low-frequency  (VLF)  and  low-frequency  (LF) 
signals  between  airborne  terminals  is  much  more  complicated  than  for 
ground-based  links,  where  only  transverse-magnetic  (TM)  waveguide  modes 
need  be  considered  for  long-range  communications.  An  antenna  at  an 
altitude  of  several  kilometers  can  radiate,  or  receive,  horizontally 
polarized  transverse-electric  (TE)  modes  with  reasonable  efficiency. 
Moreover,  for  high-speed  aircraft,  trailing-wire  antennas  have  a nearly 
horizontal  orientation,  and  are  therefore  inherently  better  coupled  to 
horizontally  polarized  signals  than  to  vertically  polarized  ones. 

Thus,  TE  modes  in  the  earth-ionosphere  waveguide  cannot  be  neglected 
for  elevated  antennas,  and  the  received  signal  will  be  a complicated 
combination  of  TM  and  TE  modes. 

This  report  analyzes  the  effects  of  antenna  elevation  and  inclina- 
tion on  the  structure  of  air-to-air  VLF/LF  communication  signals. 
Propagation  broadside  to  the  horizontal  component  of  trailing-wire 
transmitting  antennas  is  considered,  and  results  are  presented  for 

ambient  daytime  and  moderately  disturbed  ionospheric  conditions,  and 

-3  -5 

ground  conductivities  of  10  and  10  mhos/m.  For  these  ionospheric 
environments,  t le  effects  of  the  geomagnetic  field  can  be  (and  are) 
neglected.  Numerical  calculations  are  made  of  excitation  factors, 
height-gain  factors  and  attenuation  rates  of  the  first  few  TM  and  TE 
modes,  and  detailed  plots  are  given  for  each  field  component  versus 
distance  for  antenna  altitudes  between  0 and  40  kft.  The  numerical 
results  are  given  in  Sec.  II,  and  the  equations  used  are  given  in  the 
appendix. 
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II.  ELECTRIC  DIPOLE  FIELDS 


The  fields  of  an  electric  dipole  in  the  earth-ionosphere  waveguide 
can  be  expressed  as  summations  of  elementary  waveguide  modes — so-called 
TM  and  TE  modes.  For  daytime  or  disturbed  conditions,  when  geomagnetic 
mode-coupling  can  be  neglected,  the  vertical  electric  dipole  (VED) 
excites  only  TM  modes.  Each  TM  mode,  and  hence  the  total  VED  field, 

it  ' 

contains  only  three  components;  vvz,  the  radial  (.yr  vertical)  electric 

V V 

field,  E^,  the  transverse  magnetic  field,  H , and  the  longitudinal 

V 

electric  field,  EQ.  The  horizontal  electric  dipole  (HED)  excites  both 
TM  and  TE  modes  and,  for  arbitrary  orientations,  excites  all  six  field 
components.  However,  for  propagation  perpendicular  to  the  dipole  axis — 

the  situation  analyzed  in  this  report — only  the  radial  (vertical) 

H H 

magnetic  field  H , the  transverse  electric  field,  e",  and  the  longitu- 

r tp 

dinal  magnetic  field,  H^,  are  excited.  Although  these  are  the  same 
three  field  components  that  are  present  in  a pure  TE  mode,  the  HED 
broadside  fields  are,  in  fact,  composed  of  complicated  combinations  of 
both  TE  and  TM  modes.  This  absence  of  pure  modal  excitation  has  led 
to  the  nomenclature  "quasi-modes"  for  the  HED  fields. 

The  rather  complex  expressions  for  the  fields  are  given  in  the 
appendix,  as  are  the  equations  and  numerical  procedures  used  to  obtain 
the  parameters  characterizing  the  various  modes.  The  contribution  of 
each  mode  to  the  total  field  is  proportional  to  the  product  of  four 
factors;  viz , the  excitation  factor,  the  transmitter  height-gain 

* 

Conventional  spherical  coordinates  are  used,  as  are  MKS  units. 

The  transmitter  is  at  0-0,  the  waves  propagate  in  the  0 direction,  the 
path  is  at  an  azimuthal  angle,  <p,  and  the  HED  axis  is  in  the  $»0 
direction. 
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factor,  the  receiver  height-gain  factor,  and  the  propagation  factor. 

Numerical  results  for  each  of  these  factors  for  several  modes,  as  well 

as  for  the  total  VED  and  HED  fields,  are  given  below.  Earth  conductivi- 
-3  _5 

ties,  o^,  of  10  and  10  mho6/m  are  considered.  Most  of  the  results 

are  for  a nominal  daytime  ionosphere,  approximated  by  an  exponential 

conductivity  height-gradient  (B“0.5,  H - 70  km,  in  Wait's  1970  nomen- 
* 

clature;  see  appendix).  For  comparison,  some  results  are  given  for  a 
"moderate"  spread-debris  nuclear  environment  (W  * 2x10  ^ in  the  nomen- 
clature  of  Field  and  Dove , 1975,  who  give  ionization  and  collision- 
frequency  height-profiles). 

EXCITATION  FACTORS 

The  strength  with  which  a given  mode  is  excited  by  a ground-based 
transmitter  is  accounted  for  by  the  excitation  factor,  A.  Usually,  the 
VLF  excitation  factor  is  normalized  to  unity  for  TM-mode  excitation  by 
a VED  in  an  idealized,  sharply  bounded  earth-ionosphere  waveguide. 

This  normalization  is  accomplished  by  factoring  out  the  "height  of  the 
ionosphere,"  which  is  well  defined  under  the  assumption  of  a sharply 
bounded  ionosphere  ( Galejs , 1972).  This  factorization  cannot  be 
made  in  our  treatment,  since,  in  our  use  of  arbitrarily  continuous 
ionospheric  height-profiles,  the  artifact  of  an  "ionospheric  height" 
never  appears  in  the  equations.  Consequently,  the  excitation  factors 

J.  R.  Wait,  Electromagnetic  Waves  in  Stratified  Media,  Pergamon 
Press,  New  York,  1970. 

** 

E.  C.  Field,  and  M.  Dore,  VLF/LF  TE-Mode  Propagation  Under 

Disturbed  Ionospheric  Conditions , Air  Force  Cambridge  Research 

Laboratories,  AFCRL-TR-7 5-0382 , July  1975. 

***  , . , 

J.  Galejs,  Terrestial  Propagation  of  Long  Electromagnetic 

Waves,  Pergamon  Press,  New  York,  1972. 
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computed  here  have  units  of  inverse  distance,  and  differ  in  magnitude 
from  those  of  other  authors  (.Wait  and  Spies,  1964;*  Galejs,  1972)  by  a 
factor  corresponding  roughly  to  the  "ionospheric  heights"  at  which 
important  reflections  occur.  (See  appendix.)  These  heights  are,  of 
course,  ill-defined  for  a realistic  diffuse  ionosphere,  and  depend  on 
frequency,  mode  number,  and  ionospheric  model.  However,  our  excitation 
factors,  for  example,  are  smaller  than  those  of  Wait  and  Galejs  by  a 
factor  of  (6x10  - 7x10  m)  for  VLF  and  nominal  daytime  conditions. 

Figure  1 shows  the  frequency  dependence  of  the  excitation-factor 
magnitudes  of  the  first  three  TM  modes  (a -1,2,3)  and  first  two  TE 

_o 

modes  (6  - 1,2)  for  ambient  daytime  conditions  and  o - 10  mhos/m. 

g 

(A  ground  dielectric  constant  of  10  is  assumed  in  all  calculations  in 

this  report.)  The  important  thing  to  note  from  Fig.  1 is  the  relative 

magnitudes  of  the  various  modes,  since  these  magnitudes  give  the 

relative  contribution  of  each  mode  to  the  total  fields  at  short 

distances  from  ground-based  electric  dipoles.  Effects  of  elevation 

and  long-path  propagation  are  accounted  for  in  the  height-gain  and 

propagation  factors,  to  be  given  later.  At  the  lower  frequencies, 

4 4 

multiplying  by  nominal  reflection  heights  of  6 x 10  to  7 x 10  m 
gives  a result  near  unity,  as  discussed  earlier. 

The  first  three  TM  modes  (denoted  by  "a")  are  excited  about 
equally  at  the  lower  VLF  frequencies,  but,  above  about  30  kHz,  the 
higher-order  TM  modes  are  much  more  effectively  excited  than  the  first 
one.  The  TE  modes  (denoted  by  "6")  are  excited  much  more  poorly  than 

* _ 

J.  R.  Wait,  and  K.  P.  Spies,  Characteristics  of  the  Earth- 

Ionosphere  Waveguide  for  VLF  Radio  Waves,  U.S.  Department  of  Commerce, 
National  Bureau  of  Standards,  Tech.  Note  300,  1964. 
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the  TM  inodes,  as  shown  by  the  fact  that  A exceeds  A„  by  four  or  five 
orders  of  magnitude.  This  behavior  corresponds  to  the  well-known 
inability  of  transmitters  on  highly  conducting  grounds  to  effectively 
radiate  TE  waves  at  low  frequencies.  As  was  the  case  for  TM  modes,  the 
second  TE  mode  is  more  effectively  excited  than  the  first. 

The  efficiency  of  TE-mode  excitation  relative  to  TM-mode  excitation 
improves  as  the  ground  conductivity  is  reduced.  Table  1 shows  the  mag- 
nitude of  the  excitation  factors  for  the  first  two  TM  and  TF.  modes  at 
20  kHz  for  ambient  daytime  conditions  and  - 10  ^ mhos/m.  The  TE- 

mode  excitation  factors  are  nearly  two  orders  of  magnitude  greater  for 
-5  -3 

o “ 10  mhos/m  than  for  a “10  mho9/m.  On  the  other  hand,  corapari- 

o 8 

son  of  Table  1 with  Fig.  1 shows  that  the  lowest  TM  mode  is  much  more 

-5  -3 

poorly  excited  for  a * 10  rahos/m  than  for  10  mhos/m. 

5 

TABLE  1.  MAGNITUDES  OF  TM-  AND  TE-MODE  EXCITATION  FA  TORS 

AT  20  kHz  FOR  AMBIENT  DAY  AND  a = 10~5  mhos/m 

g 


1st  TM 

2nd  TM 

1st  TE 

2nd  TE 

(«-l) 

(a  * 2) 

(8-1) 

(8  * 2) 

-6 

-5 

-9 

-8 

2.4x10 

9.6x10 

9.4x10 

5.2x10 

The  excitation  factor  also  depends  on  the  state  of  the  ionosphere. 
As  indicated  above  and  in  the  appendix,  the  excitation  factors  as 
defined  here  are  inversely  proportional  to  a quantity  that  becomes  the 
"height  of  the  ionosphere"  in  the  limit  of  a sharply  bounded  iono- 
sphere. For  the  diffuse  ionospheres  used  here,  the  excitation  factors 
at  the  lower  VLF  frequencies  are  roughly  proportional  to  the  inverse 
of  the  average  height  at  which  important  reflections  occur.  Thus,  one 
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would  expect  these  factors  to  become  somewhat  larger  under  disturbed 
conditions,  where  the  ionosphere  is  significantly  depressed.  Table  2 

gives  calculated  excitation  factors  for  a moderate  spread-debris  nuclear 

-9  . -3 

environment  (i.e.,  W = 2x10  ; see  Field  and  Dove,  1975 ) and  o =10 

8 

mhos/m.  Comparison  of  Table  2 with  Fig.  1 shows  that  the  excitation  of 
the  lower  TE  and  TM  modes  at  20  and  40  kHz  is  increased  by  the  distur- 
bance. Results  are  shown  only  for  the  first  modes  because  the  higher 
modes  are  so  heavily  attenuated  as  to  be  unimportant  at  useful  communi- 
cation distances. 

TABLE  2.  MAGNITUDES  OF  LOWEST  TM-  AND  TE-MODE  EXCITATION 
FACTORS  FOR  A MODERATE  SPREAD-DEBRIS  NUCLEAR 

ENVIRONMENT,  AND  o = 10~3  mhos/m 

g 


Frequency 

1st  TM 
(a*  1) 

1st  TE 
(a  = 1) 

20  kHz 

2.4xl0-5 

6.0xl0~10 

40  kHz 

2.6xl0-5 

3.5xl0'10 

HEIGHT-GAIN  FACTORS 

The  height-gain  factor  of  a waveguide  mode  accounts  for  the 
effects  of  nonzero  transmitter  and  receiver  heights  More  precisely, 
the  transmitter  height-gain  factor  is  the  ratio  of  the  mode  excitation 
for  the  actual  transmitter  height,  h,  to  the  mode  excitation  for  a 
ground-based  transmitter.  Similarly,  the  receiver  height-gain  factor 
is  the  ratio  of  the  mode  strength  at  the  receiver  height,  z,  to  the 
mode  strength  at  the  ground.  Physically,  the  TM-mode  height-gain 
factors  are  equal  to  the  height-profiles  of  the  transverse  magnetic 
fields,  H , normalized  to  unity  at  the  ground,  whereas  the  TE-mode 

(pCt 
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height-gain  factors  are  the  height-profiles  of  the  transverse  electric 

fields,  E.fl,  normalized  to  unity  at  the  ground. 

<fo 

The  manner  in  which  the  hfei6nt-gain  factors  enter  the  field  equa- 
tions is  shown  in  detail  in  the  appendix.  But,  the  salient  point  is  as 
follows:  For  elevated  terminals,  the  contribution  of  a mode  to  the 

total  field  is  proportional  to  the  product  of  the  excitation  factor, 
the  transmitter  height-gain  factor,  and  tnt  receiver  height-gain  factor. 
The  transmitter  and  receiver  height-gain  factors  of  a given  mode  are 
functionally  identical  and,  therefore,  are  equal  when  transmitter  and 
receiver  are  at  the  same  altitude.  In  this  case,  the  relative  impor- 
tance of  a mode  Is  proportional  to  the  product  of  the  excitation  factor 
and  the  square  of  the  height-gain  factor  computed  for  the  appropriate 
altitude. 

Figures  2 and  3 show  height-gain  factors  G and  G for  the  first 

a p 

three  TM  modes  (a  -1,2,3)  and  first  two  TE  modes  (B  * 1,2)  for  ambient 

_3 

daytime  conditions,  * 10  mhos/m,  and  frequencies  of  20  and  40  kHz. 

These  height-gain  factors  exhibit  the  classic  height-dependences  for 

antennas  over  a highly  conductive  ground;  viz,  the  TM-mode  height-gain 

factors  are  of  order  unity  over  most  of  the  waveguide,  except  for  'ome 

rather  sharp  nulls;  above  a few  kilometers,  the  TE-mode  height-gain 

factors  increase  sharply  to  values  over  100. 

For  elevated  antennas,  the  large  TE-mode  height-gain  factor 

mitigates  the  effects  of  the  small  excitation  factor,  and  these  modes 

can  be  excited  about  as  effectively  as  TM  modes.  Stated  mathematically, 

for  altitudes  where  the  height-gain  factors  have  achieved  nearly  their 

2 

maximum  value,  the  quantity  AG  is  of  the  same  order  of  magnitude  as 

p p 


13 


15 


Figure  3.  Height-gai 
modes  (6=1 


the  quantity  A^G^.  Numerical  values  from  Figs.  1 through  3 can  be  used 
to  confirm  this  fact,  as  can  the  relationships 


■ <*>  2 

-1 

oc 

. f G«(r)dr 

A = 
’ 6 

2 j Gg  (r)dr 

J n2(r) 

a J 

a 
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(appendix,  Eqs.  (A26)  and  (A37)),  which  guarantees  that  large  height- 

gain  factors  are  associated  with  small  excitation  factors,  and  vice 

2 

versa.  In  the  above  equations,  a is  the  earth's  radius,  and  n is  the 
ionospheric  refractive  index.  Figures  2 and  3 show  that,  even  for 
aircraft  altitudes  of  5 to  10  km,  the  TE-mode  height-gain  factors 
become  very  large,  albeit  well  below  cheir  peak  values. 

Note  that,  at  40  kHz,  the  first  IM  mode  shows  some  tendency  toward 
the  characteristics  of  the  so-called  earth-detached  or  "whispering 
gallery"  modes,  which  have  most  of  their  energy  concentrated  at  the 
base  of  the  ionosphere.  The  case  shown  is  not  particularly  striking, 
since  the  peak  of  the  a * 1 height-gain  factor  at  45  to  50  km  is  only 
about  three.  However,  at  higher  frequencies — say,  70  to  100  kHz — the 
first  TM  height-gain  factor  can  become  larger  than  100  just  below  the 
ionosphere,  thereby  exhibiting  true  "whispering  gallery"  characteris- 
tics ( Gale.js , 1972).  The  use  of  only  a few  waveguide  modes  to  charac- 
terize 70-  to  100-kHz  waves  is  invalid  for  most  conditions. 

Figure  4 shows  the  height-gain  factors  of  the  lowest  two  TM  and 
TE  modes  for  ambient  day,  = 10  mhos/m,  and  a frequency  of  20  kHz. 
Comparison  with  Fig.  2 (p.  14)  shows  that  the  shapes  of  the  TE-mode 
height-gains  are  virtually  unchanged  by  the  two-orders-of-magnitude 
reduction  in  ground  conductivity,  but  the  peak  values  of  the 
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Figure  4.  Height-gain  factor  for  lowest  two  TM  modes  (a*l,2)  and  lowest  two 


height-gains  are  reduced  by  an  order  of  magnitude.  This  behavior  is 

consistent  with  the  two-orders-of-magnitude  increase  in  excitation 

factor  discussed  above,  since  A.  tends  to  be  proportional  to  l/o  , and 

p g 

1/2 

Gg  tends  to  be  proportional  to  o ( Galeje _,  1972).  The  main  effect  of 

g 

the  change  in  ground  conductivity  on  the  TM-mode  height-gain  is  a 
relatively  minor  change  in  shape,  which  causes  the  a = 1 mode  to  have 
slight  earth-detached  properties. 

ATTENUATION  RATES 

Modal  attenuation  rates  account  for  the  reduction — over  and  above 

cylindrical  spreading — in  strength  suffered  by  a mode  in  propagating 

between  transmitter  and  receiver.  Figure  5 shows  calculated  attenuation 

rates  versus  frequency  for  the  first  three  TM  modes  and  the  first  two  TE 

-3 

modes  for  ambient  daytime  conditions  and  a «=  10  mhos/m.  The  higher- 
order  modes  are  more  heavily  attenuated  than  the  lower-order  ones, 
which  allows  one  to  neglect  the  higher  terms  in  the  mode  sums  at  large 
distances.  At  40  kHz,  however,  the  attenuation  rates  for  the  TE  modes 
differ  by  only  about  2 dB/Mm.  Given  that  the  higher  modes  are  more 
effectively  excited  than  are  lower  ones,  it  is  doubtful  that  the  mode 
sum  can  be  truncated  throughout  the  upper  LF  band  under  ambient  condi- 
tions. In  this  case,  ray  tracing  is  a better  approach  to  the  propaga- 
tion analysis.  For  moderately  disturbed  conditions,  the  higher  modes 
are  much  more  heavily  attenuated  than  the  first  modes,  and  the  mode 
sum  can  be  used  well  into  the  LF  regime. 

Figure  5 shows  that  the  lowest  TE  mode  is  slightly  less  attenuated 
than  the  lowest  TM  mode.  This  behavior  is  in  contrast  to  results  given 
by  Field  and  Dore  {1975)  for  a perfectly  conducting  ground,  and  an 
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ambient  daytime  model  ionosphere  that  differs  slightly  from  the  one 
used  here.  Those  results  indicated  that  the  first  TM-mode  attenuation 
rate  was  between  2 and  3 dB/Mm,  which  was  1 or  2 dB/Mm  less  than  the 
lowest  TE-mode  attenuation  rate.  The  decrease  in  assumed  ground  conduc- 

-3 

tivity  from  00  to  10  mhos/m  causes  the  first  TM-mode  attenuation  to 

essentially  double,  while  the  TE-mode  attenuation  is  relatively  unaltered. 

In  addition  to  the  results  given  in  Fig.  5,  some  attenuation  rates 

have  been  calculated  for  ambient  day  with  o = 10  ^ mhos/m,  and  a 

8 

_9  _3 

moderate  (W  ■ 2x10  ) spread-debris  nuclear  environment  with  o - 10 

g 

mhos/m.  The  results  are  shown  in  Tables  3 and  4.  Comparison  of  Table  3 
with  Fig.  5 shows  that  the  TE-mode  attenuation  rate  is  relatively  insen- 
sitive to  ground  conductivity,  being  only  slightly  larger  for  o * 10  ^ 

_3 

mhos/m  than  for  a “10  mhos/m.  On  the  other  hand,  the  TM-mode 
g 

attenuation  rates  become  quite  large  for  o *■  10  ^ mhos/m.  This 

* 

behavior  agrees  with  the  results  of  Pappert  (1970) , and  occurs  because 
the  TM  reflection  coefficient  has  a quasi-Brewster  angle  effect  at  the 
ground,  whereas  the  TE  coefficient  does  not.  Comparison  of  Table  4 with 
Fig.  5 shows  that  both  TM-  and  TE-mode  attenuation  rates  are  drastically 
increased  by  the  disturbance,  with  the  TE  modes  becoming  much  more 
highly  attenuated  than  the  TM  modes.  This  behavior  is  consistent  with 
that  computed  by  Field  and  Dore  (1975)  for  infinite  ground  conductivity. 

In  summary,  TM-mode  attenuation  rates  are  more  sensitive  to  changes 
in  ground  conductivity  than  are  TE-mode  attenuation  rates,  but  TE  modes 
are  relatively  more  sensitive  to  ionospheric  disturbances. 


w .. 

R.  A.  Pappert,  Effects  of  Elevation  and  Ground  Conductivity  on 
Horizontal  Dipole  Excitation  of  the  Earth-Ionosphere  Waveguide,"  RADIO 
Science,  Vol.  5,  No.  3,  March  1970,  pp.  579-590. 
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TABLE  3.  TM-  AND  TE-MODE  ATTENUATION  RATES  AT  20  kHz 

FOR  AMBIENT  DAY  AND  a - 10-5  mhos/m 

8 


1st  TM 

2nd  TM 

1st  TE 

2nd  TE 

(o-l) 

(a -2) 

(6-1) 

(6-2) 

8.6  dB/Mm 

87  dB/Mm 

3.4  dB/Mm 

10.9  dB/Mm 

TABLE  4. 

LOWEST  TM-  AND 

TE-MODE  ATTENUATION 

RATES  FOR 

A MODERATE  C PREAD-DEBRI S NUCLEAR  ENVIRONMENT 

AND  a ” 10" 3 mhos/m 
g 


Frequency 

1st  TM 
(o-l) 

1st  TE 
(a  - 1) 

20  kHz 

11.2  dB/Mm 

24.9  dB/Mm 

40  kHz 

13.2  dB/Mm 

19.0  dB/Mm 

FIELD  STRENGTH  VERSUS  DISTANCE 


The  main  ingredients  used  in  forming  the  mode  sums — the  excitation, 
height-gain,  and  attenuation  factors — have  been  given  above.  Of  course, 
the  phase  of  these  quantities,  as  well  as  the  magnitudes,  have  been  used 
in  our  calculations.  The  formulas  used  to  calculate  the  VED  and  HED 
field  components  are  given  in  the  appendix. 

Figure  6 gives  the  electric  fie]d  components  versus  distance  for 

_3 

ambient  day,  o - 10  mhos/m,  and  a frequency  of  20  kHz.  The  super- 
g 

scripts  V and  H denote  VED  and  HED,  respectively.  The  VED  fields  are 
Independent  of  azimuth,  whereas  the  HED  fields  are  calculated  for 
broadside  ($  - tt/2)  propagation.  The  first  three  TM  modes  are  included 
in  the  VED  mode  summation,  whereas  the  first  three  TM  modes  and  the 
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first  two  TE  inodes  are  Included  in  the  HED  mode  summation.  Additional 
modes  would  have  little  effect  on  the  results  shown  in  Fig.  6. 

All  electric  field  strengths  in  this  report  have  been  normalized 
so  that 


KZ  /X  = 1 volt 
o 

where  I is  the  transmitting  antenna  current,  £ is  the  effective  antenna 
length,  X is  the  free-space  wavelength,  and  Zq  is  the  free-space  impe- 
dance. To  obtain  actual  field  strengths  in  volts/meter,  the  curves 
given  in  this  report  should  be  multiplied  by  the  actual  value  of  I£Zo/X. 

As  can  be  inferred  from  Figs.  2 through  4 (pp.  14,  15,  and  17), 
the  lowest  TM-mode  fields  and,  hence,  the  VED  fields,  do  not  depend 

strongly  on  height  for  typical  aircraft  altitudes  of,  say,  15  km  or  less. 
V V 

Thus,  Ef  and  E^  are  shown  in  Fig.  6 (and  in  subsequent  figures)  for  the 
case  where  both  transmitting  and  receiving  antenna  are  on  the  ground. 

The  fields  for  realistic  airborne-antenna  heights  differ  only  slightly 
from  those  shown.  The  HED  fields,  on  the  other  hand,  contain  TE  modes, 
which  depend  very  strongly  on  antenna  altitude.  HED  field  strengths 
are  therefore  shown  for  four  cases;  viz,  when  both  transmitter  and 
receiver  are  at  altitudes  of  0,  12.5,  25,  and  40  kft. 

The  VED  and  HED  fields  shown  in  Fig.  6 (ana  in  subsequent  figures) 
are  directly  comparable  only  when  both  transmitting  and  receiving 
antennas  have  equal  horizontal  and  vertical  components;  i.e.,  when 
tralling-wire  antennas  have  inclination  angles  of  45°.  If  trailing- 

wire  transmitting  and  receiving  antennas  inclined  10°  to  the  horizontal 

V 2 2 -2 

are  used,  E must  be  reduced  by  sin  10° /cos  10°  : 3 x 10  relative 
r 


H V H 

to  E , and  EQ  must  be  reduced  by  sin  10° /cos  10°  * 0.17  relative  to  E . 

»>  w $ * 

in  order  to  compare  received  signal  strengths. 

As  would  be  expected,  the  fields  of  the  ground-based  HED  are  orders 
of  magnitude  below  those  of  the  VED,  and  are  therefore  unusable.  The 

comparison  is  much  more  favorable  for  elevated  antennas,  however,  where 

V H 

exceeds  E^  by  a factor  of  5 or  10  at  h = 25  kft,  and  a factor  of  only 

_2 

2 or  3 at  h » 40  kft.  After  the  factor  of  3 x 10  adjustment  is  made 
for  10°  inclinations  of  transmitting  and  receiving  antennas,  the  HED 
signal  is  between  20  and  25  dB  stronger  than  the  VED  signal  for  an 
antenna  elevation  of  40  kft. 

Figure  7 is  a companion  to  Fig.  6,  and  shows  the  VED  and  HED  mag- 
netic intensities  for  antenna  altitudes  of  40  kft.  The  normalized 
intensities  H = Z^H  are  used  so  that  the  magnetic  and  electric  fields 
are  directly  comparable.  In  fact,  comparison  of  Figs.  6 and  7 shows 

that  the  results  are  essentially  identical,  provided  the  substitutions 

V V H Jl 

E -*•  H and  E -*•  c are  made.  Thus,  subsequent  results  need  be  shown 
r <j>  <$>  r 

only  for  the  electric  fields.  The  fact  that  E^  = ^ shows  that  the 

<f>  r 

HED  broadside  fields  are  dominated  by  the  TE  modes  in  the  mode  sum. 

If  loop  receivers  are  used  to  sense  the  magnetic  field,  only  the 
inclination  of  the  transmitting  antenna  need  be  corrected  for;  i.e., 
for  magnetic  reception  and  10°  wire  inclination,  the  VED  fields  need 
be  reduced  by  only  sin  10° /cos  10°  = 0.17  relative  to  the  HED  fields. 


Thus,  at  most  distances,  the  HED  and  VED  magnetic  fields  of  an  airborne 
transmitter  differ  by  a factor  of  less  than  3 for  a transmitting 
antenna  inclined  at  45°,  and  are  of  comparable  magnitude  for  an  inclina- 
tion of  10°.  See  Fig.  7. 
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Distance  (km) 


Figure  7 


Magnetic  field  strength  vs.  distance  for  elevated  VED  and 
HED  transmitters:  20  kHz,  ambient  day,  a * 10"3  mhos/m. 
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Figure  8 is  analogous  to  Fig.  6 (p.  22),  but  applies  for  a fre- 
quency of  AO  kHz  rather  than  20  kHz.  As  indicated  above,  the  HED 
results  are  not  as  accurate  as  for  20  kHz,  because  the  higher  modes  are 
not  much  more  heavily  attenuated  than  the  lower-order  ones,  and  modes 
beyond  a = 3,  8 * 2 should  probably  be  included.  Nonetheless,  the 
results  in  Fig.  8 should  be  sufficiently  accurate  for  general  conclu- 


H V 

sions.  For  antenna  elevations  of  40  kft  and  most  distances,  E.  and  E 

♦ r 


are  comparable.  Thus,  after  adjustments  are  made  for  a 10°  trailing- 
wire  antenna,  the  HED  signal  exceeds  the  VED  signal  by  about  15  to  30  dB, 
depending  on  whether  magnetic  or  electric  receivers  are  used.  Note  from 
Figs.  1 and  5 (pp.  10  and  19)  that,  at  40  kHz,  the  second  TE  mode  is 
excited  about  20  dB  more  strongly  than  the  first,  but  is  attenuated  only 

about  2 dB/Mm  more  quickly.  Thus,  out  to  about  10  Mm,  E"  as  shown  in 

<t> 

Fig.  8 is  composed  mainly  of  the  second  TE  mode. 

When  the  fields  are  dominated  by  a single  waveguide  mode,  the 
shapes  of  the  field  height-profiles  are  independent  of  distance  from 
the  transmitter,  and  are  essentially  the  height-gain  factors — suitably 
attenuated — of  the  dominant  mode.  For  multimode  propagation,  the 
situation  is  much  more  complicated  because  the  field  height-profiles 
depend  on  the  height-gain  factors  of  two  or  more  modes.  The  shape  of 
these  profiles  depends  on  transmission  range,  because  the  relative 
contribution  of  each  mode  varies  with  distance. 


To  illustrate  this  effect.  Fig.  9 shows  height-profiles  of  the 


„H 


transverse  HED  electric  field,  E , at  four  distances  from  a transmitter 


15  kft  above  the  ground.  The  curves  apply  to  ambient  day,  20  kHz,  and 


,-3 


o - 10  mhos/m.  Recall  that  a total  of  five  modes — a - 1,2,3  and 
8 
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Normalized  electric  field  magnitude 
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Distance  (km) 


Figure  8.  Electric  field  strength  vs.  distance  for  ground-based  and 
elevated  VED  and  HED  transmitters:  40  kHz,  ambient  day, 

a • 10-3  mhos/m. 


B * 1,2 — are  used  to  calculate  E^.  Yet  a comparison  of  Fig.  9 with 
Fig.  2 (p.  14)  shows  clearly  that  the  mode  structure  evolves  from  the 
second  TE  mode  at  short  distances  to  the  first  TE  mode  at  large 
distances.  Specifically,  the  profile  at  D ■=  700  km  is  similar  to  the 
B * 2 height-gain  factor;  the  profile  at  D > 5000  km  is  similar  to  the 
6 • 1 height-gain  factor;  and  the  profile  at  D ~ 1500  km  is  a mixture 
between  B * 1 and  B * 2.  A similar  evolution  with  distance  from  higher 
to  lower  modes  also  occurs  for  the  VED.  The  fact  that  the  first  TE 
mode  becomes  dominant  at  20  kHz  beyond  a few  Mm  indicates  that  a suffi- 
cient number  of  modes  has  been  retained  in  the  mode  sum.  Apparently, 

_3 

this  is  not  the  case  at  40  kHz,  ambient  day,  * 10  mhos/m,  where 
the  second  TE  mode  dominates  the  first  out  to  D ~ 10,000  km. 

Figure  10  shows  the  HED  and  VED  electric  field  strengths  for 
ground-based  and  elevated  antennas,  20  kHz,  ambient  day,  and  = 10  ^ 
mhos/m,  which  is  nominally  representative  of  propagation  over  ice. 

The  main  differences  between  these  results  and  those  given  in  Fig.  6 

_3 

(p.  22)  for  0^  * 10  mhos/m  are:  a)  Fig.  10  shows  generally  lower 

field  strengths  at  large  distances  than  does  Fig.  6,  b)  Fig.  10  shows 
V H 

larger  values  of  Eg  and  E^  at  h = 0,  which  occur  because  of  the  poorer 

—5  H 

ground  reflection  at  o * 10  mhos/m,  and  c)  E generally  compares 

8 T 

V 

more  favorably  with  E^  at  large  distances  in  Fig.  10.  Note,  however, 
that  the  long-range  results  of  Fig.  10  are  artificial,  because  no 
regions  of  10  ^ mhos/m  ground  conductivity  thousands  of  miles  in 
extent  exist  on  earth.  The  proper  calculation  would  be  a mixed 
seawater-ice  path,  which  is  beyond  the  scope  of  this  report.  The 
results  of  Fig.  10  are  meaningful  to  distances  of  about  2000  km,  which 
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Figure  9.  Profiles  of  electric  field  strength  at  several  distances  from  an 
HED  transmitter  15  kft  above  the  ground:  20  kHz,  ambient  day, 

o„  = 10"3  mhos/m. 


Distance  (km) 


Figure  10 


. Electric  field  strength  vs.  distance  for  ground-based  and 
elevated  VED  and  HED  transmitters:  20  kHz,  ambient  day, 

o = io~5  mhos/m. 


is  the  length  of  Greenland.  The  excitation,  height-gain,  and  attenua- 
tion factors  given  above  for  o =*  10  ^ mhos/m  are,  of  course,  applicable 

S 

at  any  point  having  this  ground  conductivity. 

Figureo  11  and  12  show  HEI)  and  VED  electric  field  strengths  versus 

_3 

distance  for  20  and  40  kHz,  o * 10  tnhos/m,  and  a moderate  spread- 

s 

debris  nuclear  environment.  Unlike  ambient  conditions,  only  a single 

TE  and  TM  mode  need  be  retained  in  the  mode  sum,  since  the  higher  modes 

are  very  highly  attenuated.  Thus,  mode  interference  patterns  are 

absent.  The  HED  fields  do  not  compare  favorably  with  the  VED  fields 

for  this  disturbed  environment,  and,  even  for  nearly  horizontal  electric- 

dipole  transmitters  and  receivers,  the  HED  signal  would  be  smaller  than 

the  VED  signal  beyond  a few  thousand  km.  This  behavior  occurs  because 

of  the  high  attenuation  of  TE  modes  in  disturbed  environments. 

For  convenience  in  comparing  performance  for  elevated  antennas  in 

V H 

ambient  and  disturbed  environments,  E and  F.  have  been  replotted  on 

r <p 
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Fig.  13  for  our  ambient  and  spread-debris  models,  20  kHz,  a * 10 

o 

mhos/m,  and  h = 25  kft.  The  environmental  degradation  at  the  larger 
distances  is  apparent. 
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Normalized  electric  field  magnitude 


Figure  11.  Electric  field  strength  vs.  distance  for  ground-based  and 
elevated  VED  and  HED  transmitters:  20  kHz,  moderate  spread 

debris  environment,  a_  = 10“ 3 mhos/m. 
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Figure  12 


. Electric  field  strength  vs.  distance  for  ground-based  and 
elevated  VED  and  HED  transmitters:  40  kHz,  moderate  spread 

debris  environment,  a = 10"3  mhos/m. 


Normalized  electric  field  magnitude 
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Figure  13.  Electric  field  strength  vs.  distance  for  elevated  VED  and 

HED  transmitters:  20  kHz,  ambient  and  disturbed  Ionospheres, 

a_  « 10" 3 mhos/m. 
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III.  CONCLUSIONS 

For  antennas  near  the  ground,  the  field  components  associated  with 
TM  modes  are  dominant,  even  for  electric-dipole  antennas  inclined  as 
little  as  10°  with  respect  to  the  horizontal.  For  a moderate  spread- 
debris  nuclear  environment,  the  vertically  polarized  components  again 
dominate  the  horizontally  polarized  components  for  all  reasonable 
airborne-antenna  inclinations,  antenna  elevations,  and  ranges  beyond 
about  1500  km.  This  dominance  occurs  because  of  the  high  attenuation 
suffered  by  TE  modes  when  ionospheric  reflection-heights  are  severely 
depressed.  For  ambient  conditions — and  presumably  for  weakly  disturbed 
conditions — the  horizontally  polarized  components  associated  with  TE- 
mode  propagation  are  dominant  for  antenna  elevations  greater  than 
25  kft,  and  antenna  inclinations  smaller  than  about  20°,  measured  from 
the  horizontal.  For  the  frequencies  considered,  only  slight  "whisper- 
ing gallery"  modal  characteristics  are  found. 
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APPENDIX 


MATHEMATICAL  SUMMARY 


This  appendix  summarizes  the  equations  used  to  obtain  the  numerical 
results  given  in  the  main  text.  These  equations  look  very  similar  to 
those  given  by  Galejs  (1972).  However,  this  resemblance  is  in  appearance 
only,  since  we  have  fully  accounted  for  earth  curvature  and  vertical 
inhomogeneities  in  the  ionospheric  refractive  index;  Galejs  assumed  a 
sharply  bounded  ionosphere.  Detailed  derivations  are  not  given  in  this 


appendix,  but  the  mathematical  development  follows  the  work  of  Friedman 


(1951)  fairly  closely,  and  the  method  of  Budden  (1961)  is  used  to 


numerically  solve  the  modal  equations. 

REFRACTIVE  INDEX 

For  computational  purposes,  the  ionosphere  is  completely  charac- 
terized by  specifying  the  complex  refractive  index  throughout  the  height- 
regions  that  govern  long-wave  propagation  in  the  earth-ionosphere  wave- 
guide. This  specification  requires  the  number-densities,  collision 
frequencies,  and  masses  of  each  charged  species  present.  Specifically, 
the  refractive  index,  n,  is  given  by 


n (r) 


1 - 


N 


1837e  m X^u)(w-iv  )q 
o e“™  Y Y 

Y 


(Al) 


Bernard  Friedman,  "Propagation  in  a Non-homogeneous  Atmosphere," 
Communications  on  Pure  and  Applied  Mathematics,  Vol.  4,  August  1951, 
pp.  317-350. 

K.  G.  Budden,  Radio  Waves  in  the  Ionosphere , Cambridge  University 
Press,  1961. 
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f 


where  w is  the  wave-angular  frequency,  eQ  is  the  electric  permittivity 

of  free  space,  e is  the  electron  charge,  and  m is  the  electron  mass. 

e 

The  number  density,  N^,  collision  frequency,  v^,  and  atomic  mass  number, 
q^,  of  the  y species  can  be  functions  of  altitude.  The  effects  of  the 
geomagnetic  field  have  not  been  included  in  Eq.  (Al) , since,  for  the 
daytime  or  disturbed  environments  considered  in  this  report,  its  effects 
on  mode  structure  are  negligible. 

The  analysis  described  below  is  quite  general,  and  may  be  used  for 
arbitrary  height-dependences  of  and  v^.  For  VLF  propagation  in 
ambient  daytime,  however,  only  the  electron  terms  need  be  retained  in 
Eq.  (Al).  Moreover,  the  electron  density,  Ng,  and  collision  frequency, 

ve>  may  be  accurately  represented  by  the  following  analytic  forms  over 

the  height  range  of  interest: 

Ng  - Np  exp [b (z-H) ] (A2) 

and 

Ve  “ Vo  exP [~a  Cz-lO  ] , (A3) 

where  z and  H are  in  kilometers.  Since  w <<  vg  at  the  frequencies  and 

altitudes  of  interest,  it  follows  from  Eq.  (Al)  that  the  imaginary  part 
of  the  squared  refractive  index  is  proportional  to  Ne/vg  and,  hence,  has 
an  inverse  scale  height,  B •»  a + b.  Good  approximations  to  actual  day- 
time  conditions  are  obtained  by  using  Nq  * 3.93  x 10  el/m  , vq  “ 

5 x 10^  sec  \ 6 * 0.5  km  \ and  H - 70  km.  In  the  notation  of  Wait 
(1970),  these  parameters  correspond  to  a so-called  "B“0.5,  H-70"  iono- 
sphere. For  the  disturbed  ionosphere,  actual  electron  and  ion 
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number-density  and  collision  frequency  profiles,  taken  from  Field  and 
Dore  (1975),  are  used. 

MODAL  EQUATIONS 

The  fields  of  vertical  and  horizontal  electric  dipoles  are  composed 
of  summations  of  elementary  modes — so-called  TM  and  TE  modes — that 
propagate  in  the  earth-ionosphere  waveguide.  The  most  difficult  task 
in  computing  the  fields  is  finding  the  eigenvalues  of  the  dominant 
modes.  Once  these  eigenvalues  are  known,  it  is  straightforward-  llbelt 
somewhat  tedious — to  find  excitation  and  height-gain  factors,  and  to 
perform  the  mode  summations  to  obtain  the  fields. 

A spherical-polar  coordinate  system  is  used;  the  source  is  located 
at  9 ■ 0,  and  the  waves  travel  in  the  ©-direction.  Azimuthal  (#)  sym- 
metry may  be  assumed  in  determining  the  mode  properties,  because  the 
effects  of  the  geomagnetic  field  on  mode  structure  are  small  for  the 
environments  considered.  Of  course,  the  ^-dependence  of  the  total 
fields  must  be  considered  when  mode  excitation  is  calculated  for  the 
horizontal  electric  dipole  (HED) . Throughout  the  remainder  of  this 
appendix,  the  subscripts  a and  0 will  denote  quantities  associated  with 
TM  and  TE  modes,  respectively.  Unless  otherwise  noted,  MKS  units  are 
used. 

TM  Modes 

For  TM  modes,  the  electric  and  magnetic  fields  can  be  written 

S.  ' [StEra(r'9)  + ^ <«> 

and 

2a  ' %H*a<r-e)  el“t  ' <A5) 
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For  computational  purposes,  it  is  convenient  to  define  the  normalized 


wave  admittance. 


and  the  related  quantity 


01/2  Yj 

H, 

jig 

E.. 

0a 


(A  -1) 

W = 

n ' (A  +1) 
a 


where  is  the  magnetic  permittivity  of  free  space.  By  substituting 
Eqs.  (A4)  and  (A5)  into  Maxwell's  equations  and  applying  the  constitu- 
tive relations,  it  follows,  after  some  manipulation,  that  is  governed 


by  the  following  equation: 


dW  I . a2 (1-C2) 

a k ) 2,,  „ ,2  , a , 1 

dr  = 2i  )n  (1-V  ~ 1 ~ 2 2 + 2 2 2 

/ r n 4k  r n 


1 .) 

■ C1Wa>  ( ' 

-1  ) 


where  k *=  w/c,  c is  the  vacuum  speed  of  light,  a is  the  earth's  radius, 
2 

and  n is  given  by  Eq.  (Al).  In  Eq.  (A8) , the  quantity  C can  be  inter- 
preted as  the  complex  cosine  of  the  incidence  angle  of  the  wave  at  the 
ground  (i.e. , at  r=a)  . Solution  of  Eq.  (A8)  will  satisfy  the  proper 

boundary  conditions  for  only  discrete  values,  C^,  of  C.  Hence,  is 

til 

called  the  eigenvalue  of  the  a TM  mode.  The  boundary  condition  on  W 

at  the  earth's  surface  determines  C and  is  therefore  called  the  modal 

a 

equation.  This  modal  equation  is 


W (C  ,r«a) 
a a 


n2  - Wn2  - (1-C2) 
K J g a 


n2  + Wn2  - (1-C2) 
g ' g a 
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where 


io 
_ £ 

we 

o 


(A10) 


and  e and  o are  the  conductivity  and  relative  dielectric  constant  of 
8 8 

the  ground,  respectively.  In  Eq.  (A9) , the  sign  of  the  radical  is  chosen 
to  correspond  to  a downgoing  wave  in  the  ground. 

2 

Data  for  N^,  , and  are  sufficient  to  determine  n as  a function 

2 

of  height.  Once  n is  determined,  Eqs.  (A8)  and  (A9)  form  a closed  set 
for  W and  C^,  and  are  solved  by  straightforward  iteration.  Each  itera- 
tion requires  the  numerical  integration  of  Eq.  (A8) , which  is  started  at 
a great  height  where  a purely  upgoing  wave  is  assumed  as  an  initial 
condition.  Thus,  since  the  WKB  (eikonal)  solution  may  be  used  at  great 
heights,  the  assumed  value  of  is  given  in  terms  of  the  refractive 
index,  nQ,  at  the  starting  height,  rQ,  by 


W(a,r=ro) 


n2  - Jn2-(1-C2) 
o T o a 

Vn2-(1-C2) 

lo  a 


(All) 


n2  + 
o 


where  the  sign  of  the  radical  is  chosen  to  correspond  to  an  upgoing 

wave.  In  practice,  care  is  taken  to  choose  a starting  height  so  large 

that  is  insensitive  to  its  precise  value. 

Once  C and,  hence,  W (r,C  ) have  been  determined,  it  is  a simple 
a a a r 

matter  to  calculate  all  electromagnetic  parameters  associated  with  a 
given  mode.  First,  it  is  convenient  to  define 


S = 
a 


(A12) 
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Then, 


aside  from  a geometric  spreading  term. 


H,  can  be  written 


H.  (r-a)  « exp(-ikS  D)  , (A13) 

9a  a 

where  D is  the  path  length  along  the  earth's  surface.  It  follows  from 
Eq.  (A13)  that  attenuation  rate,  n,  and  phase  velocity,  V,  are  given  by 

n - 8.7  x 106k  Im  S , dB/Mm  (A14) 


and 


V = 1 

c Re  S 


(A15) 


where  the  subscript  a has  been  suppressed.  Equations  (A13)  through 
(A15)  are  valid  beyond  a few  e-folding  distances  from  the  source  and 
not  too  near  the  antipode. 

TE  Modes 

The  TE  mode  fields  are 


h - %Vr-e)  eiwt  (A16) 

and 

Sf>  - [8rHr6(r'e)  * VWr’W]  el“'  ' <A1?) 


Define 


and 


(A18) 
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t 'hence 


6 (Ae+cg) 


k 2 2 n2  i 

_ c2a^)2-B_  X — ^~rr-2  <*-V 

/ C r n 4k  r n 


The  starting  value  for  the  numerical  integration  of  Eq.  (A20)  is 


VVr'ro} 


K-^-4 

*4 


where  the  sign  of  the  radical  is  chosen  to  correspond  to  an  upgoing 
wave  at  r . The  eigenvalue  of  the  Bth  TE  mode  is  found  from 


yvr=a) 


Vv  (1-ce>  - ce 
Vv  (i-s2) + 4 


where  the  sign  of  the  radical  is  chosen  to  give  a downgoing  wave  in  the 
earth.  The  formulas  for  attenuation  rate  and  phase  velocity  are  iden- 
tical with  Eqs.  ( A 1 2 ) through  (A15),  if  C'  is  substituted  for  C . 

fs  « 

VERT  I CAT.  ELECTRIC  DIPOLE 

”e  consider  a vertical  electric  dipole  (VED)  transmitter  located 
at  a distance,  b,  from  the  earth's  center;  h * b - a,  then,  is  the 
altitude  of  the  transmitter  above  the  ground.  The  magnetic  and 
electric  fields  of  this  transmitter  are  composed  solely  of  elementary 


TM  waveguide  modes,  and  are  given  by  the  following  equations: 


HV  - 1Ua 


JZ  JHZHY' 

l DX  ysin  D/a 


-ikS  D-Itt/4 


* '■’*s1/2^ 


b b JDX  ^slnD/a^Vo  2...  Ga(b)Ga(r^  * AT/m 

n (d; 


(A23) 
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jz  jum  v* 

™DX  ^sln  D/a  a*/ 
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1/2  e 


-ikS  D-iir/4 


IsinD/aljVa  . , % 2,.  , “a  ""''a 
a (b) 


G„(b)G~(r)  » v/m 


(A24) 


FV  ^ _ 1IZo*a  JT  J D/a  A c^/2  £ 

r 2.  v,  yDX  lain  p/a  / ^ a a ..... 

n (r)br  n (b) 


-ikS  D-iir/4 


■= G (b)G  (r)  , v/m 

Z „ . a a 


(A25) 


where  I is  the  antenna  current,  i,  is  the  effective  antenna  length,  Zq 
is  the  free-space  impedance,  and  X is  the  free-space  wavelength.  The 
TM-mode  excitation  factors  are  given  by 


A 

a 


G2  (r)dr 
a 

n2(r) 


(A26) 


and  the  height-gain  factors  are  given  by 


Ga(r) 


a 

7 exp 


_ik  f pV..’.)  dr. 

ik  I A (r’)  dr 
**  a 


(A27) 


The  admittance,  A^,  is  calculated  from  Eqs.  (A7)  and  (A8) , 
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Equations  (A23)  through  (A27)  comprise  a generalization  of  those 

given  by  Galejs  (1972,  p.  97),  which  were  derived  assuming  a uniform, 

sharply  bounded  ionosphere  at  a height  Hq  above  the  ground.  The  main 

advantage  of  the  treatment  given  here  is  the  allowance  for  an  arbitrary 

ionospheric  refractive-index  height-profile  in  the  solution  of  the 

2 

modal  equation,  as  discussed  above.  Moreover,  the  terms  Aa(r),  n (r) , 

2 

and  n (b)  in  Eqs.  (A23)  through  (A27)  account  for  the  fact  that,  at 
great  transmitter  or  receiver  heights,  the  fields  are  affected  by  the 
refractive  index’s  not  necessarily  being  unity  in  the  vicinity  of  the 
antennas.  Note  that  no  "ionospheric  height"  appears  in  Eqs.  (A23) 
through  (A29),  because,  for  arbitrary  ionosphere  height-profiles,  an 
"ionospheric  height"  is  an  artifact  that  need  not  be  defined  in  the 
treatment  used  here. 

It  is  easy  to  show  that  Eqs.  (A23)  through  (A27)  reduce  to  the 
standard  form  for  the  often-used  approximate  model  wherein  the  earth  is 
assumed  flat;  the  earth  is  a perfect  reflector  of  infinite  conductivity; 
and  the  ionosphere  is  a sharply  bounded,  perfect  reflector  of  zero  con- 
ductivity located  at  a height,  hQ,  above  the  ground.  In  this  approxima- 
tion, the  refractive  index  at  the  transmitter  or  receiver  must  be  unity, 

2 

since  the  entire  waveguide  cavity  is  a pure  vacuum.  Thus,  n (r)  - 
2 

n (b)  * 1 in  Eqs.  (A23)  through  (A27).  Also,  the  assumption  of  a flat 
earth  implies  b/a  « r/a  *=  1.  It  is  well  known  that  the  height-gain 
functions  are  of  the  form  cos[  (2a-l)  (tt /2 ) z /ho ] in  this  approximation, 
whence  Eq.  (A26)  gives 

\n  ->  1/Hq  , n f 1 . (A28) 
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Equation  (A27) , which  is  a well-known  result  for  VLF  propagation 
in  idealized  sharply  bounded  waveguides,  indicates  that  our  excitation 


factor  is,  by  necessity,  normalized  differently  than  that  defined  by 

Galejs  {1972).  Galejs  factored  the  "ionospheric  height"  out  of  his 

excitation  factor,  which  is  thus  dimensionless  and  of  order  unity. 

This  factorization  is  not  possible  in  our  generalized  treatment,  since 

the  artifact  of  an  "ionospheric  height"  never  enters.  The  result  is 

that  the  excitation  factors  defined  by  Eq.  (A26)  are  smaller  than  those 

4 5 

of  Galejs  by  a factor  between  5 x 10  and  10  . 

HORIZONTAL  ELECTRIC  DIPOLE 

The  fields  of  a horizontal  electric  dipole  are  much  more  complex 
than  those  of  the  VED,  since  both  TM  and  TE  modes  are  excited.  More- 
over, the  degree  of  excitation  depends  on  the  orientation  of  the  dipole 
to  the  propagation  path.  Here,  we  will  assume  a dipole  at  r ■ b,  8 ■ 0, 
with  its  axis  in  the  $ « 0 direction.  Thus,  $ is  the  angle  between  the 
dipole  orientation  and  transmission  path  direction. 

The  fields  excited  by  an  HED  are  given  by  the  following  equations: 
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where  the  propagators  for  the  TM  and  TE  modes  are  given  by 
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All  symbols  have  been  defined  earlier  in  this  appendix,  except  for  the 
TE-mode  excitation  factor  and  height-gain  function,  which  are  given  by 
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Equations  (A29)  through  (A38)  represent  generalizations  of  the  relations 
given  by  Galejs  (.1972,  p.  99)  in  that  they  account  for  arbitrary  height- 
gradients  of  the  ionospheric  height-gradients. 

For  arbitrary  values  of  $ , all  six  field  components  are  radiated 
by  the  HED.  For  broadside  orientation  (<p  * ir/2) , only  H , E , and  tiQ 

IT  y K7 

are  radiated;  whereas,  for  end-fire  orientation,  E , E_,  and  H,  are 

r o <p 

radiated.  However,  even  though  these  orientations  produce  the  same 

field  components  as  exist  in  TE  and  TM  modes,  the  fields  are  not  pure 

modes.  For  example,  even  for  sin  4>  = 1 , Eqs.  (A29)  through  (A31)  show 

that  only  H is  composed  solely  of  TE  modes.  E,  and  Hn  contain  contri- 
r <p  w 

butions  from  both  TE  and  TM  modes.  Thus,  broadside  fields  from  an  HED 
are  often  called  quasi-TE  modes  and  end-fire  fields  are  called  quasi- 
TM  modes. 
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